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5-h ■ Abstract 

In the investigation of the evolution of cosmological perturbations in inflationary uni- 
verse models the behavior of perturbations during the reheating stage is the most unclear 
point. In particular in the early reheating phase in which a rapidly oscillating scalar field 
dominates the energy density, the behavior of perturbations is not known well because their 
evolution equation expressed in terms of the curvature perturbation becomes singular. In 
this paper it is shown that in spite of this singular behavior of the evolution equation the 
Bardeen parameter stays constant in a good accuracy during this stage for superhorizon- 
scale perturbations except for a sequence of negligibly short intervals around the zero points 
of the time derivative of the scalar field. This justifies the conventional formula relating the 
amplitudes of quantum fluctuations during inflation and those of adiabatic perturbations 
at horizon crossing in the Friedmann stage, except for possible corrections produced by the 
energy transfer from the scalar field to radiation in the late stage of reheating. It is further 
shown that outside the above sequence of time intervals the behavior of the perturbations 
coincides in a good accuracy with that for a perfect fluid system obtained from the original 
scalar field system by the WKB approximation and a spacetime averaging over a Hubble 
horizon scale. 
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§1 Introduction 

In the standard inflationary universe model the seeds of the present large scale struc- 
tures of the universe are generated from quantum fluctuations of an inflaton field during 
inflationlJ B 101. Though the amplitude and the spectrum of the seed fluctuations are eas- 
ily calculated when a model of the inflaton is given, it is not a simple task to determined 
the amplitude and the spectrum of the corresponding large scale perturbations when they 
reenter the Hubble horizon because one must trace the evolution of the perturbations from 
the inflationary stage to a recent time. 

Of course it is not so hard to trace the evolution of perturbations during the inflationary 
stage and the simple radiation-dominated Friedmann stage after inflation. It is because 
the Bardeen parameter Z, which is defined as a linear combination of the curvature pertur- 
bation <£> and the velocity perturbation V, is conserved with a good accuracy for growing 
modes on superhorizon scales if the equation of state of matter is regular and the entropy 
perturbation is sma Hence, if the Bardeen parameter is also conserved in the re- 

heating stage which connects the inflationary and the post-Friedmann stage, it is possible 
to determine the amplitude and the spectrum of the perturbations at horizon crossing just 
by the values of the Bardeen parameter for the seed fluctuations at the inflationary stage. 
For example, it is the case in the simplification that the reheating is instantaneous or no 
large entropy perturbation is produced by reheating processes^! 1 . Actually in most of the 
work on the cosmological perturbations in the inflationary scenario this simplification is 
assumed. However, reheating is a rather complicated stage in realistic models and entropy 
perturbations are not necessarily negligible. In such situations the conservation of the 
Bardeen parameter for them is not well established. 

For example, in the early stage of reheating, the energy density of the universe is 
dominated by a rapidly oscillating inflaton field. In this stage, if one neglects the tiny 
radiation contribution, the equation of state of the inflaton shows an anomalous behavior 
and the Bardeen parameter diverges periodically when the time derivative of the inflation 
field vanishes. Accordingly the evolution equation for the Bardeen parameter, or the closely 
related evolution equation for the curvature perturbation becomes singular periodically^ 1 . 
Though this divergence is replaced by periodic large peaks if the contribution of radiation 
is taken into account and the height of the peaks decreases as the inflaton decays, the 
coefficients of the evolution equations for Z and $ still show bad behavior. Hence it is 
not far from obvious whether the values of the Bardeen parameter before and after the 
reheating stage coincide with each other or not. 

The main purpose of this paper is to show that in spite of this peculiar behavior of the 
Bardeen parameter its value in the stage dominated by a rapidly oscillating inflaton field 
actually coincides with that in the preceding inflationary stage except for a sequence of 
negligibly short intervals around the zero points of the time derivative of the inflaton field, 
and that the behavior of the curvature perturbation coincides with that for a perfect fluid 
system obtained by the WKB approximation and a spacetime averaging from the original 
scalar field system. The main idea is to use the gauge-invariant perturbation variable 
for the scalar field perturbations introduced by Mukhanov in another context^ 1 , whose 
evolution equation is completely regular and exactly solvable in the long- wavelength limit. 

The paper is organized as follows. First in the next section the exact solutions to the 
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evolution equation for the Mukhanov variable in the long-wavelength limit is given and 
their behavior is analyzed from the point of view of the parametric resonance. Further 
by converting the evolution equation into an integral equation with an iterative structure 
with respect to the wavenumber, the behavior of the solutions for general long- wavelength 
perturbations is described. The singular nature of the transformation between the curva- 
ture perturbation variable $ and the Mukhanov variable in the long-wavelength limit is 
also pointed out. Then in §3 the behavior of the curvature perturbation and the Bardeen 
parameter is analyzed in detail by solving the iterative integral equation for the Mukhanov 
variable, first for the massive free scalar field and then for a scalar field with a potential 
of a generic power-law type. In §4 the results obtained by this analysis is compared with 
the behavior of perturbations for a perfect fluid system obtained by the WKB approxima- 
tion and a spacetime averaging from the original scalar field system via a relativistic virial 
theorem. Section 5 is devoted to summary and discussion. 

Throughout the paper the natural units c = % = 1 are adopted and SttG is denoted as 
k 2 . Further the notations for the perturbation variables adopted in the article ] are used 
and their definitions are sometimes omitted except for those newly defined in this paper 



§2 Basic perturbation equations and their general anal- 
ysis 

In this section we derive regular gauge-invariant equations for scalar perturbations of a real 
scalar field and the gravitational field, and give iterative expressions for their solutions. 

In this paper we assume that the scalar field couples the gravitational field minimally 
and its Lagrangian density is given by 

= ~ {Q^d^dA + 2U{<j>)) . (2.1) 

We will not specify the potential U (0) in this section. 

Throughout the paper we only consider the case the unperturbed background spacetime 
is described by a spatially flat Robertson- Walker metric 

ds 2 = -dt 2 + a 2 dx 2 . (2.2) 
Hence the unperturbed background of the scalar field follows the equation 

+ 3^0 + ^ = 0, (2.3) 
and the scale factor a is determined by 

H 2 :=(l) 2 = ^(<P 2 + 2U). (2.4) 
From these equations the time derivative of H is given by 

H = -p 2 . (2.5) 

Throughout the paper it is understood that all the perturbation variables represent Fourier 
expansion coefficients of perturbations of the corresponding quantities around this back- 
ground. For simplicity we omit the wave number suffix k for them. 
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2.1 (X, $)-system 

The most important geometrical quantity describing scalar perturbations is given by the 
gauge-invariant variable $ defined by 

$:=-Jl-—a g , (2.6) 

where TZ and a g are the .perturbations of the three curvature and the shear of each constant 
time slice, respectivelyo Q 1 . This quantity represents the curvature perturbation of the 
Newtonian slice a g = 0. Therefore a natural gauge-invariant variable constructed from the 
perturbation 5<f) of the scalar field is 

X:=8<t>-^XT B . (2.7) 

From the components (0,0), (0,j) and the traceless part of (j, k) of the perturbed 
Einstein equations, we obtain the following evolution equations for X and 

X0-0X + 2 $ = 4^$, (2-8) 
«r or 

k 2 

$ + = -_0X (2.9) 

These equations correspond to the Hamiltonian and the momentum constraints, respec- 
tively. The equations obtained from the other components do not give new equations 
independent from these. 

On the other hand the perturbed field equation for gives© 1 

X + 3HX + + X = -40$ + 217+$. (2.10) 

This equation is also automatically satisfied under Eqs.( |2.8|) -( |2.9| ) due to the Bianchi iden- 
tity if we require smoothness of X and $. 

By eliminating $ from Eq.( p.l0|) with the help of Eqs.( |2.8| )-( |2.9| ), we obtain a closed 
2nd-order equation for X: 

X + 3HX+(^ + U H - 2k 2 A X = ^.^Xl^ m - 4>X). (2.11) 

Similarly by eliminating X, we obtain an equation for 

1 ' wV+f-^-^l* = 0- ( 2 - 12 ) 



\a0 2 / \a 2 4> 2 2 

Both of these equations are singular. Hence, though the independent degree of freedom 
is two in the current system, it is difficult to analyze the behavior of perturbations by a 
2nd-order differential equation in terms of X and $. In particular this singular nature 
make it impossible to construct solutions iteratively with respect to k, as we will see later. 
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In spite of this difficulty the above set of equations has one good point. It is that we 

can easily solve them in the long- wavelength limit k — 0. To see this, let us rescale $ as 

u:=^$. (2.13) 
Then for k = 0, Eqs.(|T8| ) -([2.9| ) are written as 

(?)'=-!* (2 - l4) 

K 2 - k 2 acb , 

" - 2H* U = (2 ' 15) 

By eliminating X from these equations we obtain 



W 2 / 2a 
This equation can be easily solved to give 



H u) =^-u. (2.16) 



u = A-^B [ ^dt, (2.17) 



2 J t0 H 2 

where A and B are integration constants, and to is some initial time. Hence for k = the 
exact solutions for (X, $) -system are given by 

These solutions are completely regular and will play an important role in the following 
analysis. 

With the help of these solutions and the Green function method we can transform the 
2nd-order equation for $ above to the following iterative integral equation: 

U2 ^ „ n7 Orr f U X 



<5> = AU 1 + BU 2 + 2k 2 U Y ( -^<S>dt - 2k 2 U 2 I —i-Qdt, (2.20) 

Jt a 2 d) 2 J t a 2 d 2 



where 



H 

U x := -, 2.21 

a 

H r ad) 2 2 / H r \ , 

It is clear that this integral equation yields an expression which is singular at dj = when 
one tries to solve it iteratively. Hence the system (X, $) is not appropriate for investigating 
the long-wavelength behavior of the perturbations as commented above. 
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2.2 Y-system 

The above difficulty in the (X, $)-system is resolved if we use a new gauge-invariant variable 
introduced by Mukhanov'i Q 1 . It is defined by 



Y:=X-— $. (2.23) 
H 



Since it is expressed as 



Y = 6<f>-^11, (2.24) 

it represents the perturbation amplitude of the scalar field on the flat slices 5TZ = 0. 
This quantity has also a close relation to the Bardeen parameters' defined by 

I 2 k 2 \ aH , 

C: =( 1 + 9(TT^wj*-- l/ ' (2 ' 25) 

where w — P/p, and V is a gauge-invariant quantity representing the velocity perturbation 
in the isotropic gauge. This gauge-invariant quantity represents the curvature perturbation 
of uniform- Hubble slices'! ] and is conserved approximately for superhorizon perturbations 
if 1 + w is a positive quantity of order unity. In this paper instead of this original Bardeen 
parameter we use the gauge-invariant variable defined by 

Z:=$- ^V, (2.26) 
k 

and refer it as the Bardeen parameter because it coincides with the original one for super- 
horizon perturbations under the same condition on 1 + w as above, and shares the same 
good property^ 1 . 

Further it is more suitable for the analysis of the problem we are concerned. In partic- 
ular in the present case in which matter consists only of the scalar field V is expressed in 
terms of X as0 ] 

V = — rX. (2.27) 

acj) 

Hence Z is simply written in terms of Y as 

Z = -S-Y. (2.28) 



By a simple calculation it can be shown that Eqs. (|2.8|) - (|2.9|) leads to the regular closed 
2nd-order differential equation of Y, 

Y + 3HY +(^ + U H + 3k 2 <P 2 - + 2k 2 ^U^) Y = 0. (2.29) 



This equation is equivalent to the original system ( |2.8|) and (|2.9| ) for k ^ because by the 
inverse transform from Y to $ and X given by 

K 2 a 2 H 



2 k 2 



(2.30) 



X = Y + — $ (2.31) 
H 
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we can recover the original equations for X and $. 

In the limit case k = 0, however, this inverse transform is singular. Hence if we take 
the limit by keeping $ finite, we obtain the additional constraint on Y given by 



<pY -Hi 



H 



Y = 0. 



(2.32) 



Due to this the exact solutions ( [2.18| ) and ( |2.19|) for the X and $ in the limit k = reduce 
to the single solution for Y, 



Y = B- 



H' 



(2.33) 



It is easily confirmed by a direct calculation that this satisfies Eq. Q2.29|) for k = 0. Of 
course all the solutions for k = can be constructed from this special solution with the 
help of the standard procedure as 



Y = AV\ + BV 2 ; 
Vx:-- 



H 



H 2 



V 2 := — — ^dt. 
H J to a?(t) 2 



(2.34) 
(2.35) 

(2.36) 



The new independent solution V 2 is apparently singular at <fi = 0. However, it is 
completely regular actually as is expected from the regularity of the original equation for 
Y. To be more precise, it has a finite limit at = and can be uniquely extended across 
that point smoothly. In order to find an explicit expression for this extended solution V 2 , 
let us denote the sequence of time at which = by Tj, and examine the behavior of V 2 
around t = tj. 

If we take a sequence of time tj so that Tj_i < tj < Tj, then V 2 should be expressed in 
each interval tj_i < t < Tj cLS 



V 2 = AaVi + Bj— [ -^dt 



(2.37) 



where A = and B = 1. In order to find the recurrence relations for Aj and Bj, let us 
expand at t = Tj cLS 



(2.38) 



where the suffix j represents the value at t — Tj. Inserting this expansion into Eq. (|2.3| ) and 
Eq.( [2.4| ), we obtain 



-U, 



4>h 



-3H 



H = Ha 



l + OUt-Tj 



a = aj 



l + H^t-T^ + O^t-Tj) 2 ) 



(2.39) 
(2.40) 

(2.41) 
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Hence the integrand of V 2 is expressed as 

H 2 



<j) 2 a 3 4> 2 a 3 



(t - Tj y 



+ 0(1) 



(2.42) 



From this it follows that in the interval Tj-i < t < tj the integration in Eq. fl2.37f ) is 
written as 



H 2 

u cb 2 a 3 



-dt 



H 2 



o 2 a'> o]<v)U ■ r,) 



H 2 

dt + 



1 



1 



(2.43) 



and in the interval Tj < t < r j+1 as 



H 2 

t j+ i <fi 2 a 3 



-dt 



3+1 



H 2 



<j) 2 a 3 4> 2 a 3 (t - tj) 



j H j 
dt+ 3 



jtf \ 7 j 1 T j - tj+i. 



(2.44) 



Since the integrands in the right-hand sides of these equations are regular at t = Tj, V 2 has 
the finite limits at t = Tj ± given by 



V 2 {Tj - 0) 



V a ( Ti + 0) = -S 



J. .V 
(pjClj 

Hj 

3+ 1 T~3 m 



'<f>j a j 



(2.45) 
(2.46) 



Hence from the continuity of V 2 at t — Tj we obtain 

Bj = B = 1. 
Similarly V 2 also has the finite limits given by 



(2.47) 



Hj <p 3 a 3 \2 



3 tf,- 1 



T j -tj. 



Hj u 



H 2 



<J) 2 a 3 0|a] (* - T S ) 



dt, 



V 2 (T j + 0)=Aj +1 ^- + 



H 



Hj 0j (i i 



3 V2 3 



T j ~ tj+i , 



H 



H 2 



<J) 2 a 3 $af (t - Tj) 



(2.48) 



dt. 



(2.49) 



Hence the continuity of V 2 gives the recurrence relation 



A j+1 - Aj 



H 2 



<j) 2 a 3 <f) 2 a 3 (t - Tj) 



dt- % 



1 1 

+ 



djO] \Tj ~tj tj +l -Tj / 



(2.50) 



In the next section we will show that we can take tj so that it is near the zero points of 



and Aj = 0. 



7 



This argument shows that V2 takes a definite value independent of to at t = Tj, 

V ^ = ~ ( H ^( v (2 - 51) 

which changes its sign alternatively with j. We will show in the next section that it 
approximately gives the amplitude of V2 at each interval and its absolute value decreases 
monotonically with j. Hence V2 is a decaying mode. On the other hand from Eq. ( |2.4|) the 
amplitude of V\ stays nearly constant. This behavior should be contrasted with that of 
4> whose amplitude decreases monotonically in the rapidly oscillating phase, and can be 
understood result of the parametric resonance. 

For example, for the case U = m 2 <j) 2 /2, the equation for Y is written as 



i 3 / 2 Y)- + 



2 4m 
m H — — sin 2m(t - t c 



i 3 ' 2 Y) ~ 0, 



(2.52) 



under the condition (k 2 /a 2 )/mH <C 1 and mt ^> l(see Eq. (|3.2| )). This equation satisfies 
the condition for the parametric resonance to occur. Since it can be put into the integral 
equation 



fl 3/2y = Acosm{t - t ) + B sin m(t - t 

4 



+ cos m(t — to) [ dt sin m(t — to)- sin2m(t — to)a 3 ^ 2 Y 
Jto t 

— sin m(t ~ to) / dtcosmit — t )- sin2m(t — to)a 3 ^ 2 Y, (2.53) 

Jtn t 



the amplitude of the solution corresponding to A = 1 and B = is multiplied every each 
period by 

r to+At 2 At 

dt- sin 2 2m{t - to) ~ 1 + — , (2.54) 



1 + 



to 



where At is the period of oscillation. Hence after a long time a 3 / 2 Y behaves as 

f — t 

a 3 ' 2 Y ~ e Jt o 1 cosm(t — to) = — cos m(t — to)- 

to 



(2.55) 



Similarly for the solution corresponding to A = and B = 1 the amplitude is multiplied 
by 1 — At ft every each period and behaves as 



— f — ^0 
a 3 ' 2 Y ~ e J 'o 1 cosm(t — t ) = — cos m(t — to). 

It is easy to see that this behavior coincides with that of V\ and V2 approximately. 
On the other hand Eq.( [2.3|) is written as 



(2.56) 



,3/2, 



+ 



ur ■■■()[ i 



(a 3 ^) = 0. 



(2.57) 



In this case, since dt/(mt 2 ) converges to a finite quantity of order 1 / (mto), no parametric 
resonance occurs for mto ~ 1- 



S 



Here note that the amplitude of Y increases in proportion to a power of time by the effect 
of parametric resonance in contrast to the Mathieu function in which the amplitude grows 
exponentially. This feature is common to the parametric resonance of coupled oscillating 
fields on expanding universes whose scale factor increases in proportion to a power of time 
as the above argument shows(cf. Ref.JTT] ]). 



2.3 Iterative solutions 



When the wave number k is small, we can convert the differential equation ( 2.29 ) into an 



iterative integral equation with the help of the solutions V\ and V 2 for k = and the Green 
function method: 



Y = AVi + BV 2 + k 2 Vi f aV 2 Ydt - k 2 V 2 f aV x Ydt 

= AVx + BV 2 + k 2 Go(aY). (2.58) 
From this equation and Eq. ( |2.30D $ is expressed as 

* = T^ w(yi > Y) > (2 ' 59) 

where W is the Wronskian of Eq. ( |2.29| ) given by 

W{u,v) := a 3 {uv -vu). (2.60) 

Let Y\ and Y 2 be the solutions to the above integral equation for (A,B) = (1,0) and 
(A,B) = (0,1), respectively, and $1 and $2 be the corresponding solutions for $. Then 
since 

w(y 1 ,Y 2 ) = w(y 1 ,v 2 ) = i, (2.61) 

from Eq. ( 12.59) ) $1 and $2 are written as 

(2.62) 
(2.63) 

We will show in the next section that Y\ ~ V\ and Y 2 ~ V 2 , and the second term on the 
right-hand side of Eq. (|2.63|) is much smaller than unity for superhorizon scale perturbations. 





K 2 H 








/ aVxYxdt, 








K 2 H 




$ 2 = 


1-k 2 [ a 


2 ak 2 ' 




\ Jt 



From this it follows that 



$1 ~ const, (2.64) 



Hence in spite of the singularity of the transformation (|2.30| ) and (|2.31|) the solutions ( |2.18 



and (|2.19j) for the (X, <3>)-system with k = yield a good approximation for the behavior 
of superhorizon scale perturbations if we give appropriate ^-dependences to the constants 
A and B. 
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§3 Behavior of the curvature perturbation and the 
Bardeen parameter 

In this section we solve the integral equation for Y given in the previous section and analyze 
the behavior of the curvature perturbation and the Bardeen parameter for superhorizon 
perturbations in the oscillatory phase of the scalar field. First we look at the case U is 
proportional to (f) 2 in detail and then extend the analysis to the generic case U oc \<fr\ n . 



3.1 U = m 2 (j) 2 /2 case 

In the case in which the potential of the scalar field is given by 

1 



U 



2/2 

-m <p , 



(3.1) 



Eq. 



is rewritten as 



,3/2, 



+ 



m 2 + 



H + -H 2 

2 



( a 3 / 2 0) = 0. 



From Eqs. (|2.4j) and (|2.5|) the solution to this equation is estimated as 

a 3 / 2 = a sin m(t - t ) + —sin mt [ (h + -H 2 ] cosmt a 3/2 (f>dt 

2m Jt \ 2 J 



2m 



t 



— — cos mt I I H + -H z ) sin mt a 3/2 <j)dt 



a 



sinm(t — 1 ) ~ 



1 



1 



where 



Am \to t 

_ H 

m 



cosm(t — t ) + 0(e 



Here note that in the oscillatory phase of (p, e <C 1. Similarly <fr is estimated as 



1/15 

cos m(t — t ) H I ) sinm(t — t ) + 0(e 

Am \t t 



c?l 2 <\) = ma 

By inserting these expressions into Eq. ([2.4j) and solving it we obtain 

1 



H 

„3 



2 

3f 
3 



2mt 

K 2 a 2 {mtf [l + 0(e 2 ) 



sin2m(t-t ) + 0(e 2 ) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 
(3.7) 



Next we determine the explicit behavior of V\ and V2 with the help of these equations. 
Since we need an estimate for a 3 <j) 2 /H 2 which is correct up to the order (t — Tj) at each 
interval tj < t < tj+i, we first estimate its second derivative and then integrate it. 
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Let us denote the sequence of the zero points of <fi by Tj and take another sequence of 
time, tj, as in the previous section. Then from Eq. ([3.5|) the quantity £j defined by 



satisfies the equation 



(31 



cos6 



^ ArrtTj AmtQ 



sm 



( i + 0(ef) 



where ej = efa). This gives the estimate 



sin^ = (-iy + e 



Hence in the interval tj < t < tj+i, (a?/ 2 cf)/H) is estimated as 



(3.9) 



(3.10) 



"' ,3/21 



1 + 2k 1 K 



9H 2 k 4 



{-if-cmfCj 



where 



H 4 m? 4 m? 2 m?H 2 
9 

mrj sin a; + x sin a; + 7 (cos 3a; — cos a;) + O(xe ) 



x := m(t — Tj), 



and Cj is a constant depending on j such that 



a -, = 1 + e 



0; ' 



(3.11) 
(3.12) 

(3.13) 



which is related to the value of (a 3,/2 0/if) at t = Tj by 



fa) 



'a 3 / 2 m 2 N 



fa) = (-1) J+1 ^W' 



(3.14) 



The integration of this equation gives the following estimate for the integrand in Eq.( [2.36| ): 



H 



SIM 



ttit. 



+ f(x) + O(x 2 e ) 



where 



f(x) :— x sin 2x. 



(3.15) 



(3.16) 



Inserting these expressions into the definitions of V\ and V2, and using the recurrence 
relation (|2.50|) , we obtain 



V 2 



'~ iy+1 ^ [mT J + f{x) + ° (eo)] Sinx ' 



3mcra 3 / 2 



1 + 



rriTi 



mTj 



2x 



mTj 



cosx 



(3.17) 
(3.18) 



11 



From Eq. ( |3.9| ) sin x and cos x are written in terms of t — to as 

5 



(— l) 3 sin x = cos m{t — t ) 
(— iy cosx = sin m(t — to) + 

Hence Eqs.( p.l7|) - (|3.18|) are written as 
V6 



Amt Amt 
5 1 



— ymm(t-t ) + o(el), (3.19) 



Amt Amt 



cosm(t-to) + 0(e^J . (3.20) 



Vi 



Vo 



K 



cosm(t — 1 ) + 
AV2 1 



1 



1 



Amto mt 



sin m(t — to) + 
1 



sin m(t — t ) H sin 3m(t — t ) + O (ef) 

Amt v ' 

(3.21 

3 1 



Amt 



2mt Amto 
cos3m(t — t ) + O^e^ 



cos m(t — 1 ) 



These expressions show that in the oscillatory phase we can take tj so that 

sinra^ - t ) = O(e ) , 
V 2 {t 3 ) = 0. 



(3.22) 



(3.23) 
(3.24) 



They also show that V 2 is approximately in proportion to 0/t, and takes its local minima 
and maxima at around t = Tj. Further the amplitudes of V\ and V 2 are estimated as 



Vl = O 

v 2 = o 



K,m 3 a 2 t 2 ) 



(3.25) 
(3.26) 



Hence V 2 is a decaying mode. 

With the help of these estimations we can easily estimate the behavior of Y in the case 
in which k is small but nonzero by solving the integral equation fl2.58j ) iteratively. First 
integrals appearing in the first-order iteration are calculated as 



to 



aV{dt 



2 ^ - la ~ (/,//„)' 



aVxV 2 dt 



AOa^m 
7 



1 - 



tn\ 5 / 3 



t 



+ 



to 



lQaom 2 
2 



l + O(6 )- (^) 4/3 {l-^sin^m(t-t 



6mt 
2 



sin2m(t-t ) + 0(eg)(3 T 27) 

(3.28) 



to 



aVndt - 



1 + e 



to \ 7/3 
t 



7 



21agm 4 a 2 t 

Hence in the second order in k, Y\ and Y 2 are given by 



6mt 



sin2m(t — tc 



(3.29) 



Y l {t)= V^t) 
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3amt e / to 



10a 3 / 2 11 V t 



1/3 



tn\ 5 / 3 



l-lf) +O(eo) 



sinm(t-t ) + 0(V) , (3.30) 



, 4 e 



o(i 4 ) , 



63 Zq ma 2 (mt ) 2 



Vi(t) 



O(e ) 



(3.31) 



where I denotes the ratio of the perturbation wavelength and the Hubble horizon radius, 

aH 



k 



(3.32) 



In order to obtain the estimate in full order in k, we use an inequality obtained from 
Eq. fl2.58|) instead of the iterative power series. Let us consider the integral equation 



y = v + k 2 G o y. 

If we introduce \y\ m (t) and \v \ m (t) defined by 

\y\ m (t) := sup \y(t')\, 
t <t'<t 

\v\ m (t) := sup \v(t')\, 
to<t'<t 

the integral equation yields the inequality 

\y(t)\ < \v\ m {t) + k 2 \V 2 \ f a\Vi\dt + \Vi\ f a\V 2 \dt 

J to J to 

Since the integrals in the right-hand side of this equations satisfy the inequalities 

9(7 



\y\m(t). 



a\V!\dt< 1/2 
*o i(Ja m 



(mto^imtf 3 [1 + O(eo)] 
2 1 



J a\V 2 \dt<^j 2 -^--[l + O(e )] 



(3.33) 

(3.34) 
(3.35) 

(3.36) 

(3.37) 
(3.38) 



from Eqs. (^2lD -(^22D, we obtain 



8 Co 

\y(t)\ < \v\ m (t) + -4\y\ m (t). 

6 < 



(3.39) 



Since |u| m (t) and \y\ m (t) are nondecreasing functions, we can replace \y(t)\ by \y\ m if) in 
this inequality. Hence we obtain the estimate 



\y(t)\ < \y\ m (t) < 



\v\m(t) 

1 - 3e / " ; 



(3.40) 



Now let us apply this inequality to Y\ and Y 2 . First we rewrite the integral equation 
for Y u 

Yi = Vi + k 2 G o Y\, (3.41) 
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as 



Y 1 -V 1 = k 2 GoV 1 + k 2 Go(Y 1 -V 1 ), (3.42) 

and apply the inequality (|3.40| ) by putting y = Y\ — V\ and v = k 2 G o V\. Then since the 
first-order estimate Eq.( |3.3U| ) gives 



k 2 \GoVA<^[l + 0{e»)\ 



we obtain the inequality 



1^1-^1 < 



V / 66 / 2 [l + O(e )] 



Similarly from 



k 2 \Go V 2 \ < 



5k 1 — 3eo^ 2 
2 e 1 



[1 + O(e )] 



21 a^ 2 m 2 at 
the deviation of Y 2 from V 2 is bounded as 

2 1 eo /- 2 [l + 0(e )] 



\Y 2 - V 2 \ < 



21 a^ 2 m 2 at 1 — 3e / 2 



(3.43) 

(3.44) 
(3.45) 

(3.46) 



With the help of these estimate we can easily determine the behavior of the curvature 
perturbation and the Bardeen parameter. First the curvature perturbations $i and $2 
corresponding to Y\ and Y 2 are given by 



$ 5 



I - [^ 5/3 + O(e ) + of| 




k 2 H 



2 ak 2 



(3.47) 
(3.48) 



Hence $2 decreases in proportion to t~ 5 ^ 3 while $1 approaches -3/5 for t >i . 

On the other hand the Bardeen parameter is simply related to Y as Z = —HY/<j>. Hence 
it oscillatorilly diverges. However, its value is well-controlled and stays almost constant 
except for very small intervals around t = Tj in the situation corresponding to the reheating 
phase after inflation. 

The inflationary phase, during which is very small, terminates when e = H/m be- 
comes smaller than unity and the scalar field begins to oscillate. Since Z stays constant 
during the inflationary staged 1 , and since \&/Z\ is of order <fi 2 /p and is much smaller than 
unity at the end of the slow-rolling phase^ 1 , the behavior of the perturbations in the oscil- 
latory phase after the inflation is well described by the solution Y with the initial condition 
$(' r o) — and Z(t ) = Z Q , where Z is the value of Z during the inflationary stage. From 
Eq. ( |2.28 ) and Eq. ( |2.30|) this initial condition is expressed in terms of Y as 



Y(t ) = 0, Y(t ) = -jtZ , 



(3.49) 



where the suffix denotes the values at t — tq. 
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If we write Y as the linear combination of the fundamental solutions Y\ and Y 2 in the 
form 

Y = AY 1 + BY 2} (3.50) 
the constant A is expressed in terms of Z as 



A = -W(Y 2 ,Y) = ^Z Y 2 (r ), 

-no 

since W{Y X , Y 2 ) = W(V U V 2 ) = 1. From Eq. (^51|) and Eq. (|3~46l) Y 2 (r ) is given by 



(3.51) 



r 2 (r ) = -^ + O| 



3/2 72 
, a am l o . 



(3.52) 



Hence we obtain 



Similarly from 



A = - 



1 + 



£0 



Zn. 



(3.53) 



B = W(Y 1 ,Y) 



Hn 



Z Y 1 (t ) 



Y 1 (t ) = O[ 



leo\ 



(3.54) 
(3.55) 



we obtain 



5 = 



' cr 2 m 2 ' 



Hence from Eqs.( |3.47D -( |3.48| ) the behavior of $ is determined as 



(3.56) 



$(£) = A$i(f)+B$ 2 (t) 



-Z 



1 - 



to \ 5/3 
£ 



O(e ) + O^ 



fo 
*8, 



(3.57) 



Thus $ stays constant for £ ^> £ . 

On the other hand, in the interval tj < t < £j+i, the Bardeen parameter Z = —HY/<p 
is given by 



Z(t) = Z 



1 + o| U01 



1 



eo (to \ 1/3 



sinm(£ — Tj) Zq V £ 



(3.58) 



Hence at least except for the small intervals m|£ — Tj\ < 1/Iq, Z coincides with the constant 
Z with the accuracy O(eo/Zo), which is of order 10~ 20 for perturbations relevant to the 
present large scale structures. 
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3.2 U = X\(J)\ n /n case 

Let us extend the analysis in the previous subsection to a more general case in which the 
potential is given by 

U=^\(l>\ n , (3-59) 

where n is a positive constant equal to or greater than 2. 

First note that from the equations for the background quantities 

+ 3#0 + A|0| n /0 = 0, (3.60) 

the condition for to oscillate, \H<j)\ <C A|0| n_1 is written as 

K\(f>\ < 1. (3.62) 

Under this condition, from 

= 00 + 3#00+ A|0|" 

= (m> + ' - (l - ^ 2 2 ) 2 + nU(4>), (3.63) 

we obtain the cosmic virial theorem 

< U(<j>) >= - < <p 2 > [l + 0(e 2 ) + O(HT)] , (3.64) 
n l v / j 

where < Q > denote the time average of Q over the cosmic expansion time scale, T is the 
period of oscillation and 

e := K < 2 >V2 (3. 65 ) 

which is of the same order as e introduced in the previous subsection for U = m 2 2 /2. 
From this virial theorem < p > is expressed in terms of < <fi 2 > as 



" 1J r l + 0(e 2 )+0(OT)]. (3.66) 



<p>= ^r < ^ 2> 



Hence from the energy equation 

p = -3H<P 2 , (3.67) 
we obtain the following evolution equation for < p >: 

< p >■ = -^^H < p > [l + 0(e 2 ) + 0(HT)] . (3.68) 

Solving this equation, we obtain 

< p >oc a (3.69) 

aoct 2 ^, (3.70) 

<0 2 >oca"^2. (3.71) 
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To be precise, these equations hold with the accuracy 0(e). 

Using these equations we can determine the behavior of the perturbations by the same 
method as in the previous subsection. First, since the fundamental solutions V\ and V 2 for 
k = are bounded as 



K 



\v 2 \ = o 

the integrals of aVjVk are calculated as 

aV{dt 



k 2 (n + 2)(2n+ 1) 
/ a\V x V 2 \dt = / 

•/to 

f aV?dt = 

Jtn 



(at - onto) [1 + O(e )] 



a 2 H 2 



n 2 e 2 



^a 5 H 3 J 

Hence the deviation of Yj from Vj is estimated as 

1 e 

' k e 2 



1*1 


-Vi\ 


= 0( 


1*2 


-v 2 \ 


-°( 



0(K)0 



From this estimate the behavior of $j are determined as 



$1 
$2 



3n 



2(2n + 1) 



a 3 Hl 2 
3rn 
ao^o 



0(^2) 0(- 



+ 0(e Q ) + O 



k 2 H 



1 + 0, F 



2(2n + l) 
OC t 3n 



(3.72) 
(3.73) 

(3.74) 
(3.75) 
(3.76) 

(3.77) 
(3.78) 

(3.79) 
(3.80) 



Hence $2 gives a decay mode while $1 approaches a constant for t 3> to, as in the harmonic 
case. 

In the situation considered at the end of the previous subsection these estimates lead 
to the same conclusion. In fact for the solution Y satisfying the initial condition (|3.49|) , A 
and B are given by 



A = -Z 
B = 



1 + I 



Zn. 



Hence we obtain 



Z(t) 



3n 



2(2n+ 1) 



O 



£0 

A, 



at 
+ 



+ 01 



e 



+ 0(r n ) + ()(-)+o(J 



(3.81) 
(3.82) 

(3.83) 
(3.84) 
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for tj < t < tj+i. Thus Z(t) coincides with Z with the accuracy 0(y/e/l) except for tiny 
periods H\t — Tj\ < 0(y/e/l). 

§4 Comparison with the system obtained by WKB 
approximation 

In this section, after showing that when averaged over a spacetime region of the Hubble 
radius scale the behavior of the spacetime structure and the energy-momentum tensor of a 
nearly homogeneous oscillating scalar field on an expanding universe is well described by 
a perfect fluid, we compare the behavior of the curvature perturbation and the Bardeen 
parameter of that fluid system with those considered in the previous section. 

In order to find the averaged behavior of the energy-momentum tensor of a scalar field, 
let us extend the cosmic virial theorem for an exactly homogeneous scalar field derived in 
§3.2 to a nearly homogeneous case. The argument is quite similar to that there. First we 
multiply the field equation 

□0-^ = 0, (4.1) 

by y^g~(f> to get 

=v /3 W> n 0-v /3 W>^> 

= ^(v^0V» - ^Tg ((V0) 2 + 4>Ut) . (4.2) 

By integrating this equation over a spacetime region of a proper size L, we obtain 

< v^(V0) 2 > + < V=ff^ >= O (j- < I V0| >) . (4.3) 

Hence if the parameter e defined by 

i<V=?MIW|> 



L < v^(V0) 2 > 



(4.4) 



is much smaller than unity, we obtain the following relativistic virial theorem for the 
potential U = \\4>\ n /n: 

< ^~gU >= - l - < v^(V0) 2 > [1 + 0(e)]. (4.5) 

Note that e is of the same order as that in the previous section if L is taken to be 1/H. 
Applying this virial theorem to the energy-momentum tensor of the scalar field, 

T£ = V>V,0 - % Q(V0) 2 + C/) , (4.6) 

we obtain 

_ _ 77—2 

< v/=£Z? >=< ^-gg^d^dA > < v^W) 2 > [1 + 0(e)]. (4.7) 
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If the scalar field and the spacetime is nearly spatially homogeneous, we can neglect the 
correlation among g^ v and <f> up to the linear order of perturbation: 

< V^gg^dxWvt >~ 4^99^ < <9 A <R0 >, (4.8) 

< y/=jjTZ >~V=9< T£ >, (4.9) 

where we have written < g^ u > simply as g^ u . In this approximation the spacetime average 
of T*f is given by 

n — 2 

< T£ >=< V>V,0 > +-^-^ < (V0) 2 > • (4.10) 

When the scalar field oscillates rapidly enough and the spacetime is nearly flat on the 
time scale of oscillation, it can be assumed that is well approximated by the WKB form 

<j> = F(S,x), (4.11) 

where S is a rapidly oscillating phase, and VP(the partial derivative of F with respect to 
the argument x) is much smaller than the norm of the vector V0: 

/ <(VF)»> \ f <(VF)»> \ 

In this situation the correlation between (d$F) 2 and d^S is small: 

< (d s F) 2 d^Sd u S >=< (d s F) 2 >d^<S>d u <S>+0(e< (V0) 2 >) . (4.13) 

Hence if we introduce two scalar functions p and P, and a time-like unit vector by 

(p + P)Ey/„ :=< (<9 5 P) 2 > <9 M < S > d v < S >, (4.14) 

P := ^— < (V0) 2 >, (4.15) 

the above equation for < > is written as 

< T£ >~ (p + P^C/, + <^P, (4.16) 
n — 2 

P = wp; w = . (4.17) 

" n+2 K ' 

Hence the averaged dynamics of the scalar field is described by a perfect fluid. In particular 
when the geometry and the field is spatially homogeneous, the averaged energy density and 
the scale factor behave as 

poca^+^oc-^, (4.18) 

H oc a"i (1+w) oc (4.19) 
aoct^n^J. (4.20) 
This behavior coincides with those obtained in the previous section. 
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Next let us examine the behavior of perturbations of this fluid system and compare it 
with that obtained in the previous section. For scalar adiabatic perturbations for the fluid 
with P = wp, their evolution equations are expressed in terms of the curvature perturbation 
$ and the Bardeen parameter Z asS 1 



D* + ^^$ = §(1 + «,)Z, 



DZ 



where D is the differential operator 



2 w 
$ 

3l 2 l + w ' 



D 



d_ 

da 



By noting the time-dependence of I 2 and a 3 H obtained from Eq. ( |4.19|) 

1 



/ 2 



oc a 



l+3w 



a 3 H oc ai {1 - w K 



this evolution equation leads to the following 2nd-order differential equation for Z: 

a 3 H 



D(a HDZ) 



-w- 



f- 



-Z. 



By the Green function method this equation is converted into the integral equation 

B 2w f Z da 2w 1 f a 3 H ^da 



A 



Z da 

777— + 



a 3 H 3(1 - w) Jao P a 3(1 - w) a 3 H 
From this $ is expresses as 



r- 



-z—. 

a 



9(1 -w 2 ) I 2 



B 



2w 



a 3 H „da 



I 2 



Aw a 3 H 

Let us denote the solutions to this integral equation as 

$ = A® x + ^$2, 
Z = AZ X + BZ 2 . 

Then by iteration the fundamental solutions ($j, Zj) are estimated as 

2w (l 1\ Aw 1 



1 



3(l-w)(l + 3w) V^ 2 I 
3(1 + w 



3(1- w)(5 + 3w I 2 



at 



5 + 3w 

1 



a 3 H 
4 



aH \l 2 
Aw 



+ 



2w 



3(1 -w)(9w- I) I 2 3(1 -w)(l + 3w) V^ 2 i< 
Aw 1 



3(1 -w)(9w- l)%a%H 
9(1 -w 2 ) I 2 



Aw a 3 H 



, + 0| ?) +0 (g 



4.21) 
4.22) 

4.23) 

4.24) 

4.25) 

4.26) 

4.27) 



4.28) 
4.29) 



4.30) 
4.31) 



4.32) 
4.33) 
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These expressions show that for the fluid system with w 



in 



2)/(n + 2) behaves 



exactly in the same manner as those obtained in the previous section, Eqs.( [3.47l )-( 13.48| ), 
apart from a normalization constant. 

A similar conclusion can be obtained for the behavior of the Bardeen parameter. To 
see this, let us consider the solution with the initial condition 



$(t ) = 0, Z(t ) = Z . 



(4.34) 



Then, since in the present case $(t ) = and $2(^0) 7^ 0, the corresponding solution is 
proportional to ($1, Z\). Hence it is given by 



3n 



2(2n + 1) 



a^to 
at 



+ 



Z(t) = z 



0(1 



(4.35) 
(4.36) 



These expressions coincide with Eqs. (|3.85 )- (|3.84 ) in the leading order including the co- 
efficients, though the conservation of the Bardeen parameter holds better in the exact 
treatment done in the previous section than that in the WKB fluid approximation. 



§5 Summary and discussion 

In this paper we have analyzed the behavior of scalar perturbations on superhorizon scales 
in a stage dominated by a rapidly oscillating scalar field <ft and have shown that the Bardeen 
parameter stays constant in a good accuracy expect for a sequence of negligibly short 
intervals around the zero points of <fi. We have further shown that the scalar field system 
reduces to a perfect fluid system by a spacetime averaging when the WKB approximation 
is good, and that the behavior of perturbations for that perfect fluid system coincides in a 
good accuracy with that of the original system outside the above sequence of time intervals. 

In realistic situations in which radiation is produced by the decay of the scalar field 
(inflaton), the divergence at around t = r n is replaced by a sequence of narrow peaks and 
their height decreases as the energy density of radiation increases. Hence the above result 
practically establishes the conservation of the Bardeen parameter during the reheating 
phase dominated by the scalar field, provided that the energy transfer from the inflaton to 
radiation does not affect its behavior. 

Further the above result shows that we can study the evolution of superhorizon per- 
turbations during the reheating phase of inflation with a good accuracy by replacing the 
scalar field by a perfect fluid obtained by the WKB approximation and the spacetime av- 
eraging. This provides a generalization to relativistic superhorizon perturbations of the 
result obtained by Nambu and Sasaki in the Newtonian approximation^ and will signif- 
icantly simplify the analysis of the evolution of perturbations during reheating. The effect 
of the energy transfer on perturbations during reheating will be discussed utilizing this 
simplification in a next paper. 
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